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ABSTRACT 

In this paper, we explore the plausible luminosity evolution of early-type galaxies in 
different cosmological models by constructing a set of pure luminosity evolution (PLE) 
models via the choices of the star formation rate (SFR) parameters and formation 
redshift Zf of galaxies, with the observational constraints derived from the Hubble 
Space Telescope (HST) morphological number counts for elliptical and SO galaxies of 
the Medium Deep Survey (MDS) and the Hubble Deep Field (HDF). We find that the 
number counts of early-type galaxies can be explained by the pure luminosity evolution 
models, without invoking exotic scenarios such as merging or introducing an additional 
population. But the evolution should be nearly passive, with a high Zf assumed. The 
conclusion is valid in all of the three cosmological models we adopted in this paper. 
We also present the redshift distributions for three bins of observed magnitudes in 
F814w pass-band, to show at which redshift are the objects that dominate the counts 
at a given magnitude. The predictions of the redshift distribution of 22.5 < bj < 24.0 
are also presented for comparison with future data. 

Key words: galaxies: elliptical and lenticular, cd - galaxies: evolution - galaxies: 
luminosity function, mass function - cosmology: miscellaneous. 



INTRODUCTION 



One of the most basic astronomical methods is to simply 
count the number of galaxies as a function of apparent mag- 
nitude. Such exploration can be traced back to the early 
work of Hubble (Hubble 1926), and revived in the seventies 
after a nearly half century's gap (Brown & Tinsley 1974). 
Ever since then, it has been widely-used to probe the evo- 
lutionary history of galaxy populations, or to help estimate 
the parameters of cosmological geometry. It is just these 
attempts that brought about a problem called the excess of 
faint blue galaxies (FBGs), which remains one of the grand 
astronomical issues for a long time (Koo & Kron 1992, Ellis 
1997). The difficulties lie in that one cannot find a logi- 
cally simple and self-consistent way to explain the obser- 
vational data of different aspects. If one tries to reproduce 
deep blue galaxy counts using flat faint-end luminosity func- 
tions (LFs), with no evolution (nE) assumed, the observed 
number counts show an excess with respect to the expected 
predictions by a factor of ~ 4 to 5 at bj ~ 24, rising up 
by a factor of ~ 5 to 10 at bj ~ 26 (Maddox et al. 1990a; 
Guiderdoni & Rocca-Volmerange 1991), and continuously 
increasing to the faintest levels observed at bj ~ 28. The 
adoption of an open geometry can moderate, but still not 
help ameliorate the situation. The introduction of pure lu- 
minosity evolution (PLE), which, as a natural consideration, 



allows the spectra and luminosity of galaxies to change with 
redshifts, can provide a better fit to the faint galaxy number 
counts in the bj pass-band, if an open geometry (qo < 0.5) 
or a cosmological constant (A > 0) is involved (Broadhurst, 
Ellis & Shanks 1988; Colless et al. 1990; Cowie, Songaila 
& Hu 1991; Colless et al. 1993; Fukugita et al. 1990). But 
as infrared photometric data become available, one finds 
that 7f-band (2.2/^m) counts show no excess with respect 
to the no-evolution predictions up to K ~ 21 (Cowie et 
al. 1990; Gardner, Cowie & Wainscoat 1994), and such 
PLE models will overpredict the counts in Jf-band (Cowie 
1991). On the other hand, the PLE models also overpre- 
dict a high- z distribution of galaxies, which is not found in 
the results of z surveys of faint galaxies (Broadhurst et al. 
1988; Colless et al. 1990; Koo & Kron 1992 and references 
therein). To get away from such a dilemma (i.e., the opti- 
cal/infrared and the photometric/spectroscopic paradoxes), 
a number of less straitforward models concerning number 
evolution of galaxies have been proposed. One is the merger 
model which would decrease the comoving number density 
of objects while increasing their luminosities (Guiderdoni 
& Rocca-Volmerange 1991; Broadhurst, Ellis & Glazebrook 
1992; Carlberg & Chariot 1992; Kauffmann, Guiderdoni & 
White 1994). Another is to introduce an entirely new pop- 
ulation of dwarf galaxies which once existed in early times 
and have faded and/or disappeared in recent epochs (Broad- 
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hurst ct al. 1988; Cowie 1991; Cowie ct al. 1991; Babul & 
Rees 1992; Babul & Ferguson 1996). 

Another complexity we need to consider is that there 
may exist some inherent uncertainties in the present-day 
LFs, which is not well-determined by the local surveys. For 
instance, if we adopt a high normalization of the charac- 
teristic density, or consider a steep faint-end-slope, which 
may accommodate more dwarf galaxies than the flat one, 
the degree of excess will be substantially decreased (Saracco, 
Chincarini & Iovino 1996). An extensive investigation of lit- 
erature shows that there exist many discrepancies between 
different research groups (cf. King & Ellis 1985; Loveday 
et al. 1992; Marzke et al. 1994a, 1994b; Im et al. 1995; 
Roche et al. 1996) in the determination and/or adoption of 
the values of the luminosity function (LF) parameters ((f)*, 
M* , a) and the relative mix ratio between different mor- 
phological types. Some authors (De Propris et al. 1995) 
even announced that an extreme steep slope of LF faint- 
end had been observed, say a ~ -2.2, in the cores of four 
rich clusters of galaxies (Abell 2052, Abcll 2107, Abell 2199 
and Abell 2666). The uncertainties in local LFs include, 
e.g., a large local fluctuation ( Shanks 1989), significant lo- 
cal evolution (Maddox et al. 1990b), selection effects and/or 
incompleteness (Zwicky 1957; Disney 1976; Ferguson & Mc- 
Gaugh 1995), or perhaps systematic errors in local surveys 
(Metcalfe, Fong & Shanks 1995). Thus the universality of 
LFs is doubtable, while the well-determination of present- 
day LF is of great significance for understanding the galaxy 
evolution at high redshifts, and is conducive to reconciling 
the paradoxes mentioned previously. 

Great progress of observational cosmology in recent 
years has been made through the powerful ground-based 
10 — m Keck telescope, and especially the Hubble Space 
Telescope (HST). The latter, with its high resolution of 0.1" 
FWHM, can provide us with image information of great 
value in that the morphologies of different Hubble types 
can be segregated into several wide classes (cf. Driver & 
Windhorst 1995a; Driver et al 1995b; Abraham et al. 1996, 
and references therein). With morphological data, it be- 
comes possible to simplify the modelling of FBGs so that 
each morphological type can be allowed to be modelled in- 
dependently, and hence the complexity of each individual 
model can be greatly reduced (Driver & Windhorst 1995a). 

Following this line of thought, we consider first the mod- 
elling of E/SO's number counts in the Medium Deep Sur- 
vey (MDS) and the Hubble Deep Field (HDF) obtained by 
WFPC2's F814w pass-band on board HST in the present 
paper, by means of PLE models according to the latest ver- 
sion of population synthesis models of Bruzual & Chariot 
(1997, hereafter BC97), without invoking the more exotic 
scenarios mentioned above. We choose the parameters of 
the star formation rates (SFRs) and the formation redshift 
(zj) of galaxies to reproduce the photometric properties such 
as colors and spectral energy distributions (SEDs) of the 
local population. Another motivation of this paper is to 
explore up to what extent the parameters of cosmological 
models within the framework of standard cosmology can be 
constrained by our PLE models. We find that the number 
count-magnitude relation of E/S0 galaxies can be well ex- 
plained by our PLE model in any cosmological geometry we 
adopted in this paper, including a) fio=l, h=0.5 (Ho—100h 
km s _1 Mpc -1 ), b) O =0.1, h=0.5 and c) fi =0.2, A =0.8, 



h — 0.6, under the appropriate choice of the parameters of 
SFRs and z/s. This seems to be conflicting with the result 
of Drive et al. (1996). 

In Section 2 we will demonstrate the basic methods for 
modelling the number counts by pure luminosity evolution 
of galaxies and describe the cosmological models as well as 
the LF for early-type galaxies adopted in this work. Some 
details such as dust extinction are also introduced. The 
results of our models are shown in Section 3 for comparison 
with the observational data. In Section 4, we will discuss the 
influence of uncertainties in LF, dust extinction and different 
initial mass function (IMF) upon our results. A comparison 
of our work with others is also made in this section. We will 
give our summary and conclusions in Section 5. 

2 MODELLING THE NUMBER COUNTS OF 
ELLIPTICALS 

The number distribution of galaxies between the interval 
(m, z) ~ (m + dm, z + dz) is stated by the expression as 
follows: 

d 2 N{m x ,z) = ^L^L ( f ) (M x )dm x dz, (1) 
4-7T dz 

where lj is solid angle, dV is the comoving volume element 
of the redshift interval z ~ z + dz, and <f> is the present-day 
LF. The relation between apparent magnitude m and M in 
the A-pass band is given by 

m x = M x + 5/05(^| ) + corr, (2) 

where dz, is the luminosity distance measured in pc and de- 
pendent on cosmology. The term corr refers to the correc- 
tion needed to translate the galaxy observer-frame magni- 
tude into its rest-frame magnitude. In the nE model, it is 
only characterized by fe-correction, which accounts for the 
red-shifting of the spectra due to the recession of galaxies. 
The e-correction should be considered if the intrinsic galaxy 
luminosity evolution is involved, and in more realistic cases, 
some secondary effects, e.g., dust extinction, should also be 
taken into account to some extent. The magnitude-limited 
number counts (z-distribution) can be obtained by integrat- 
ing (1) over the specific magnitude range (mi, m-z). By 
integrating over redshift z the differential number counts 
(number-magnitude relation) of galaxies can be obtained. 
Since galaxies will become drastically faint beyond the red- 
shift zl at which the Lyman continuum break is shifted into 
the effective wavelength of the pass-band being considered 
(Madau 1995), the upper limit of the integral over z should 
be z up =mm(zf , zl), where Zf is the assumed z of galaxy 
formation. For B-band zl is 4.0, while for F814w it is 7.8. 
We will see the significance of this consideration in what 
follows. 

2.1 Cosmological Models 

As one of the goals of our present exploration, we would ex- 
amine the validity of our PLE models under the three choices 
of the currently popular cosmological models. The first we 
consider is Oo = l and //o=50km s _1 Mpc -1 , a flat Einstein- 
de Sitter universe favoured by the standard inflation theory. 
The second we adopted is the Friedmann-Robertson- Walker 
model, in which we adopt fi =0.1, 7io=50 km s _1 Mpc" 1 , 
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Table 1. Parameters of the cosmological models un- 
der consideration, as well as the assumed formation 
rcdshifts of galaxies. 



Scenario 


h a 


n 


Ao 




A 


0.5 


1.0 





3.5 , 5 , 8 , 10 


B 


0.5 


0.1 





2.5 , 3.5 , 5 , 10 


C 


0.6 


0.2 


0.8 


2.5 , 3.5 , 5 , 10 



a H = WOh km s" 1 Mpc" 1 ; 

b Zf refers to the redshift of galaxy formation. 

representing the low density universe. The third is A- 
dominated Friedmann-Lemaitre model for which we choose 
Ao=0.8, f2o=0.2 and a higher Hubble constant fz~o=60km 
s _1 Mpc -1 , also compatible with the inflation theory. Both 
the low density and the A-dominated universes have larger 
volumes at a given redshift. We list the parameters of these 
cosmological models in Table 1. Also listed in it are the 
galaxy formation redshifts we assumed in the three cosmo- 
logical models. 

Once a specific cosmological model is given, the re- 
lations between redshifts, look-back times and volume el- 
ements, which are needed in modelling, can be deter- 
mined completely by the standard formulae (cf. Weinberg 
1972; Guiderdoni & Rocca-Volmerange 1990; Fukugita et al. 
1990). 

2.2 Evolutionary SEDs of Galaxies 

There are standard evolutionary synthesis techniques (Bruz- 
ual & Kron 1980; Guiderdoni & Rocca-Volmerange 1987; 
Chariot & Bruzual 1991; Bruzual & Chariot 1993, hereafter 
BC93; Bruzual & Chariot 1997, hereafter BC97) to obtain 
the SEDs of galaxies. Our galaxy SEDs are computed on 
the basis of the latest galaxy isochrone synthesis spectral 
evolution library (GISSEL96, BC97). The BC97 models are 
built from a library of stellar tracks in the theoretical H- 
R diagram, covering all evolutionary stages of stars with 
multi-metallicity. The spectra of galactic stars, from near- 
UV to near-IR, extending to the far-UV by means of model 
atmospheres, are employed in the synthesis. For the present 
purpose, however, we do not model the evolution of metal- 
licity with respect to z, and consider the solar metallicity 
only. 

Colors of galaxies, as a significant photometric prop- 
erty, can provide us with valuable information about their 
present-day composition of stars, and accordingly give us 
important clues to the formation and evolution of galaxies. 
In the present work, we use several broad-band colors of 
elliptical galaxies (see Table 2) observed locally as the pre- 
liminary constraints on our selections of the model param- 
eters. The photometric systems we use here are: Johnson's 
UBVRI systems (Johnson & Morgan 1953), Couch's bj and 
r f (Couch & Newell 1980) and the Palomar 200 IR detec- 
tors' K band filter (effective wavelength, 2.2/im). We set the 
photometric zero point by the SEDs of the star of a Lyre. 
Thus we can compute the model colors to make comparison 
with the observed ones, and among them we select the SFR 
parameters to give the best fit. 

Throughout this work, we adopt the Scalo (1986) IMF 
for our models. The Scalo IMF is less rich in massive stars 



Table 2. Observed colors of local ellip- 
tical and SO galaxies for reference. 



CI 


observed 


U - B a 


0.43 


B-V a 


0.95 


B - R b 


1.83 


bj - I c 


2.39 


bj - K b 


4.16 



a Mixed from Fukugita et al. (1995); 
b Mixed from Yoshii et al. (1988); 
c Mixed from Yoshii et al. (1988) and 
Lidman et al. (1996). 



than the Salpeter (1955) IMF because of the steeper slope of 
the former at the high-mass end. By the adoption of Scalo 
IMF, the UV flux of early times can be greatly reduced so as 
to avoid a large number of galaxies being detected at high 
z, which are not observed in current deep surveys (Pozzetti, 
Bruzual & Zamorani 1996, hereafter PBZ96). 

2.3 Dust Extinction 

Simple luminosity evolution models (e.g., PBZ96) only con- 
sider the contributions to the fluxes by star populations. 
But in reality, most galaxies are full of interstellar medium 
(ISM), and hence effects of internal absorption by dust 
are also needed to be introduced into the modelling. In 
particular, the dust extinction has great influence upon z- 
distributions of galaxies (Wang 1991; Gronwall & Koo 1995; 
Campos & Shanks 1995). 

Unlike the situation in spiral galaxies, however, physi- 
cal and evolutionary relationships between the various com- 
ponents of the ISM in ellipticals are not yet well under- 
stood (Goudfrooij, 1996). To circumvent the complexity 
and uncertainty in determining the law of extinction for el- 
lipticals, we make an ad hoc assumption that Wang's pre- 
scription, which is just the simulation for spirals, could be 
extrapolated to the description of dust extinction for ellip- 
ticals, except that we choose a larger geometrical param- 
eter £=0.50, and a smaller optical depth r z= o = 0.10, in 
B-band for present-day L* ellipticals, which mimic the ge- 
ometrical feature and less dust content for local ellipticals, 
in contrast with the adoption by Wang, that £ = 0.25 and 
r z= o = 0.20 for spirals. We also assume that the optical 
depth depends on the galaxy luminosity as r cx L° z = , which 
is a simple representation of the fact that luminous galaxies 
are seen to be much dustier than those of lower luminosity. 
The extinction curve of our models is taken to be a power 
law in wavelength, namely oc A~ n , with n = 2 (Draine & 
Lee 1984), which is the same as Wang's adoption. By such 
choices, the extinctions for B-band are 0.05 mag at present- 
day, and 1.02 mag at z = 5; for F814w they are 0.02 and 
0.50 at 2 = and 5, respectively. Observationally, extinc- 
tion corrections in B band for ellipticals in the nearby 4 
clusters are 0.18mag (A2256, 2=0.0601), 0.07mag (A2029, 
2=0.077), 0.05mag (Coma, 2=0.023), and 0.04mag (A957, 
2=0.045), obtained from the NASA Extragalactic Database 
(cf. Schade, Barrientos & Lopez-Cruz 1997) . Hence, the re- 
sults of our models can be accepted when compared with the 
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Figure 1. Extinction in the B-band as a function of redshift 
for an L* galaxy. The predictions of models are shown by lines, 
with the model parameters (see Wang 1991 for denotations of 
these parameters) exhibited in the figure. 

observation. We show the relation of extinction against red- 
shift in B-band in Figure 1. It can be seen that present-day 
extinctions are not affected by the geometrical parameter (, 
but are sensitive to the local optical depth r z= o; while £ will 
play an important role at high redshifts, at which fluxes are 
greatly reduced, and hence galaxies become much fainter by 
a thick geometrical parameter than the thin one. 

Needless to say, it is only a simple treatment, since we 
do not model the dust evolution and some other details such 
as the inclination of galaxies are also neglected. But even 
so, the effect of dust extinction can greatly reduce the UV 
flux in early epochs due to the wavelength dependence on 
the extinction law. 

2.4 Luminosity Function 

The present-day luminosity function of galaxies is well rep- 
resented by the Schechter (1976) analytic form as follows 

<A( L)dL = 0*(^)V^*d(^), (3) 

where L* is the characteristic luminosity and a is the faint- 
end slope. <f>* is the characteristic density, which is the nor- 
malization related to the number of luminous galaxies per 
unit volume (cf. Ellis 1997). 

We adopt the morphology-dependent LF from the mod- 
els of Roche et al. (1996, hereafter RSMF96) for our present 
study. For early-type galaxies, the parameters of LF as- 
sumed in RSMF96 are M* B = -21.00, a = -0.70 and 
<t>* = 9.68 x 10~ 4 Mpc~ 3 . We derive the local LF of 7 F8 i4 
band by shifting M B according to the present-day color 
B — /f814 computed by our models (see the last column 



in Table 3). Note that in RSMF96, the E/S0 galaxies are 
divided into two types, namely, the cold and the hot ellipti- 
cals. In our case, however, we do not make such a distinction 
and simply incorporate them into a single type . Besides, it 
should be mentioned that the classification between early- 
type, late-type spirals and irregulars of RSMF96 LF is also 
suitable for future investigations, although we only deal with 
ellipticals in the present work. For PLE models, in contrast 
to the cases of number evolutions, the evolution of LF with 
respect to z is realized only by e-corrections ( with the dust 
extinction involved), which are functions of z. 

Considering that there exist some uncertainties in de- 
termining the local LF, as mentioned in the introduction, 
we attempt to vary the parameters of LF, in particular, a, 
to show to what extent such uncertainties can affect our 
conclusions. 

3 RESULTS 

3.1 Star Formation Rates 

As outlined previously, we adjust the model parameters 
of SFRs to fit best the observational local colors of E/S0 
galaxies. The SFR is a single-parameter function taking 
the form of exponential decrease with respect to time, i.e., 
ip(t) oc exp(—t/T e ), which is a natural result under the sim- 
ple assumption that %l>(t) is proportional to the available gas 
(Kennicutt 1983), where r e is the e-folding time characteriz- 
ing this form of SFRs. Thus in the viewpoint of modelling, 
the evolutionary photometric and spectroscopic properties 
of galaxies can be completely determined by these parame- 
ters in our simple PLE models once the cosmological models 
are specified. We summarize the values of these model pa- 
rameters in Table 3 in column 2 and 3. The model colors 
of U - B, B -V, B - R, bj - I and bj - K are listed from 
column 4 to 8, for comparison with the observational ones 
shown in Table 2. The bj — If»14 colors are also given in 
column 9 for converting the LF from bj into 7fsi4 band. 

Table 3 shows that the modelled colors in a specific 
band computed by any one of our models are much close 
to each other, the difference from model to model is very 
small. The model colors are also close to the observed ones 
listed in Table 2. In particular, the BC97 models can repro- 
duce better the bj — K color than the old version (BC93) by 
about 0.2 — 0.3mag reddening. Considering that the uncer- 
tainties in local colors of galaxies are about 0.1 — 0.2mag, 
the adopted values of the parameter r e should be accept- 
able. Furthermore, we can also find that the models can 
grossly reproduce the observed SEDs for E/S0 galaxies of 
present-day by comparing the fc-corrections computed from 
the models with the empirical one (not shown). 

Superior to the old version as it is, the BC97 models still 
need further refinement, especially in the UV range, for the 
modelled U — B colors are slightly redder than the observed 
one, which signals such improvements to be needed. 

3.2 HST F814w Number Counts of Ellipticals 

The wider spatial coverage of the MDS and the improved 
signal-to-noise of the HDF, which provide us with abundant 
image information, now allow us to model the morphologi- 
cally segregated galaxy number counts, and we pay atten- 
tion to the number counts of E/S0 galaxy in our current in- 
vestigation. The observed data of elliptical galaxy number 
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Table 3. Zf and SFR parameters listed in column 2 and 3 in the three cosmological models 
under consideration, with the modelled colors listed from column 4 to 9. 



Scenario A 



Models 


z / 


T a 

e 


U-B 


B-V 


B — R 


bj-I 


bj-K 


bj — Ifsu 


Al 


3.5 


0.8 


0.64 


0.92 


1.67 


2.27 


4.10 


2.03 


A2 


5 


0.2 


0.65 


0.92 


1.69 


2.30 


4.15 


2.05 


A3 


8 


1.0 


0.64 


0.92 


1.68 


2.29 


4.13 


2.04 


A4 


10 


1.0 


0.65 


0.92 


1.68 


2.29 


4.13 


2.04 



Scenario B 



Models 


z / 


1 e 


U-B 


B-V 


B- R 


bj-I 


bj-K 


bj — Ifsi4 


Bl 


2.5 


0.2 


0.67 


0.93 


1.70 


2.32 


4.19 


2.07 


B2 


3.5 


0.1 


0.69 


0.94 


1.72 


2.34 


4.22 


2.09 


B3 


5 





0.71 


0.95 


1.74 


2.36 


4.25 


2.11 


B4 


10 


1.0 


0.71 


0.95 


1.74 


2.36 


4.25 


2.11 



Scenario C 



Models 


z / 


T a 

e 


U-B 


B-V 


B- R 


bj-I 


bj-K 


bj — Ifsu 


CI 


2.5 


0.4 


0.70 


0.93 


1.72 


2.33 


4.20 


2.07 


C2 


3.5 


0.2 


0.70 


0.94 


1.73 


2.35 


4.23 


2.09 


C3 


5 


0.05 


0.72 


0.95 


1.74 


2.37 


4.26 


2.11 


C4 


10 


1.2 


0.71 


0.95 


1.74 


2.36 


4.25 


2.11 



a r e is the characteristic time-scale for SFRs, measured in Gyr, and r e = represents the single 
burst model, in which the evolution is pure passive. 



counts are taken from Glazebrook et al. (1995) and Driver 
ct al. (1995b) for MDS, and Abraham et al. (1996) for 
HDF. Although there may exist some errors in the morpho- 
logical classification (the scheme and details in classification 
have been discussed in the above-quoted works), from Fig- 
ure 2 we can see that the data show gross agreement within 
the overlapped regions between different research group's, 
indicating, to some extent, that the classification is reli- 
able. Another noteworthy feature is the flattening beyond 
Ifsu ~ 22.5, which is just the unique property that is dif- 
ferent from the counts of other types (spirals or irregulars) 
as well as the overall population. 

Figures 2- (a), (b) and (c) show differential number 
counts derived from the F814w bandpass on board HST, 
together with our models in the three scenarios A, B and C 
indicated by lines. Apparently, the non-evolutionary predic- 
tions cannot reproduce the data, either at faint magnitudes 
or at the bright end. The predictions with luminosity evolu- 
tion, however, can match better the number counts at bright 
end up to If814 ~ 21.0mag, no matter which scenario, and 
which Zf in each scenario is taken into account, i.e., the 
number counts at bright magnitudes are insensitive to ge- 
ometry or Zf . The effect of luminosity evolution is significant 
and should not be neglected. Furthermore, allowing for the 
fluctuation of the Glazebrook et al. sample at bright mag- 
nitudes, we believe that the normalization of the RSMF96 
LF employed in this work is reasonably well. 

Discrepancies between these evolutionary predictions 
do exist and occur beyond If814 > 21.0mag. We will discuss 
them in detail the the following. 



3.2.1 Scenario A 

All the predictions fall within the error bars of the data, 
while only the model of Zf =5 predicts better the number 
count than the others. Besides, it can grossly reproduce 
the faint-end flattening at around 7f814 ~ 23.0 — 25.0. For 
the case of Zf = 3.5, though it keeps on increasing up to 
~ 23.0mag and reproduces better the counts even than that 
of Zf = 5 model at 7fsi4 ~ 21.0 — 23.0, it will turnover 
consequently and greatly underpredict the faint-end counts 
as well as the flattening. It seems that Zf = 3.5 for a fiat 
Einstein-de Sitter world model is not proper, since it is too 
low, and hence there is not enough volume to accommodate 
sufficient objects. While if the Zf is chosen as high as 8 or 10, 
the models will predict much lower counts from 7fsi4=20.5 
to 22.5 than the others. 

3.2.2 Scenario B and C 

Models of Zf — 10 can be absolutely ruled out: the predic- 
tions at 7^814 > 24.0 overshoot the HDF data by about a 
factor of 3 to 4, indicating that volumes in the two world 
models at such a redshift are both too large; they simulta- 
neously underpredict the counts at Jf814 = 21.0 — 24.0 mag 
(Scenario B) or 21.0 - 23.0 mag (scenario C), indicating that 
luminosity evolutions in the two cases are both too low. 

The prediction seems better for Zf = 2.5 in Scenario B, 
since it can reproduce well the counts almost at every mag- 
nitude from 21.0 to 25.0 mag. It is also not bad for Scenario 
C, though the prediction at 23.0 - 24.0 mag is slightly higher 
than that of Scenario B. 
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Figure 2. The differential number counts for E/SO galaxies as a function of apparent magnitude in /f?8l4 band. Figures 2-(a), (b) 
and (c) are for Scenario A, B, and C, respectively. The sources of the observational data are exhibited in the figure and these data are 
indicated by symbols. The models are shown by lines. 
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The predictions are not as good as those of Zf = 2.5 for 
models of Zf = 5.0 in either Scenario B or C. In particular, 
the slopes 7 = dlogN(m)/dm of the models at 7fsi4 > 22.5 



mag are ~0.22 and ~0.18 for Scenario B and C, respectively, 
steeper than the flattening of the observational data. It, 
however, should be attributed not to the evolution but to the 
geometry, since the evolution is nearly passive (see Table 3). 
As can be seen, the slopes even for non-evolutionary models 
are still as high as ~0.12 and ~0.14 for Scenario B and C, 
respectively. Regardless of the steep slopes, there is neither 
sufficient nor convinced evidence to exclude such models as 
Zf — 5.0. After all, our models reproduce the counts within 
the magnitude range where there exist observational data 
(/fsi4 = 17.5 - 25.0 mag). 

The Zf — 3.5 models of Scenario B and C can be re- 
garded as interpolations between the models for Zf = 2.5 
and Zf = 5, but they are more similar to those of z/=5. 



3.3 Redshift Distribution 

In modern astronomy, redshift distributions of galaxies are 
the key measurement required to probe the local luminos- 
ity function, or possible evolution taking place in galaxies 
as well as the geometry of the universe. Such measurement 
has already reached up-to-date as faint as B ~ 24 mag by 
the new spectroscopic samples derived from the Keck Tele- 
scope (Cowie et al. 1995, 1996). It can be seen from Section 
3.2, that the conclusions are rather uncertain and obscure. 
Hence, we should incorporate redshift distributions for bins 
of observed magnitudes, to continue our investigation. Such 
exploration will be conducive to the understanding of our 
models, though such spectroscopic data has not been ob- 
tained morphologically as yet. 
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Figure 3. The z-distribution of F814w < 21.0. Panels (a), (b) and c are for Scenario A, B and C respectively. Models arc indicated 
by lines. The model predictions have been normalized in the unit of deg — 2 with redshift bin Az=0.01. 



3.3.1 z- distribution of lFsi4-band 

(a) J F814 < 21.0 

In Figure 3 we plot our predictions for the redshift dis- 
tributions of Jf814 < 21.0 mag. Panels (a), (b) and c are for 
Scenario A, B and C, respectively. We re-scale our predic- 
tions in the unit of deg -2 with redshift bin Az=0.01. We can 
see that all the models with evolution will peak at z ~ 0.5 
(hereafter, we call it the first), but it can also be seen that 
there is another peak existing at the redshift of Zf (hereafter, 
we call it the second), except for Scenario A. In particular, 



the second peak of z/ = 5.0 model in Scenario B is more than 
two times as high as the first. An extremely rapid declining 
of the SFR (for the Zf = 5.0 model in Scenario B, the time 
scale T e = 0, i.e., single burst) will consume (almost) all the 
gas in galaxies to form stars (almost) in an instantaneous 
burst, rendering the galaxies to be detectable at that high 
redshift, not bringing about the subsequent star formation, 
especially of massive stars, and resulting in the rapid fad- 
ing in luminosity of galaxies shortly after the formation of 
stars of the first generation. Therefore, the existence of the 
second peak should be completely attributed to the rapid lu- 
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Figure 4. The z-distribution of 21.0 < F814w < 23.0. Panels (a), (b) and c are for Scenario A, B and C respectively. Models are 
indicated by lines. The normalization is the same as Figure 3. 



minosity evolution. Yet the peak is very narrow, and hence 
it will not be expected to contribute too much to the counts. 

(b) 21.0< Ifsi4 < 23.0 

Figure 4 is plotted for the predictions of the redshift 
distribution limited in 21. 0< If&u < 23.0, the others are 
the same as Figure 3. For Scenario B and C, it can be seen 
that the high-z peaks of the z/ = 5.0 models are even higher 
than those in the case of If&u < 21.0, but the contributions 
to the counts are still not comparable to the first. As for 
Scenario A, it is obvious that the high-z contribution begins 



to be significant for the models whose z/=3.5 or 5. What 
forms the sharp contrast is that there are no such peaks 
in the predictions of the z/=8 or 10 models, It is just the 
absence of these peaks that leads to the lower number counts 
at around 7fsi4 ~ 20.5 to ~22.5 mag than those of z/=3.5 
or 5.0 models (see Figure 2-(a)). 

(c) 23.0< Jf814 < 25.0 

Figure 5 is for the redshift distributions predicted by 
our models within the magnitude range 23. 0< Jfsi4 < 25.0. 
All the curves of the models are similar to those in the case 
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Figure 5. The z-distribution of 23.0 < F814w < 25.0. Panels (a), (b) and c are for Scenario A, B and C respectively. Models are 
indicated by lines. The normalization is the same as Figure 3 and Figure 4. 



of 21. 0< Jfsi4 < 23.0, except that high-z contributions be- 
come dominant in such faint magnitude range, especially for 
Scenario A, as can be seen from the figure. 

In each of the three figures above-mentioned, we have 
also plotted predictions of non-evolutionary models of z/ — 
5.0 for each panel to make comparison with the evolution- 
ary predictions. The most significant difference between the 
non- and the evolutionary models can be seen from these 
figures is that no high-z peaks appear in the predictions 
of non-evolutionary models. Hence, according to the PLE 
models, and differing from nE ones, the objects at high-z 



also contribute to the predicted number counts at magni- 
tudes fainter than Jf814 ~ 20.5. Our PLE models presented 
here are characterized by these high-z peaks in the predic- 
tions of redshift distributions. 



It is worthwhile to mention that these peaks, in prin- 
ciple, are detectable through the F814w band-pass since zl 
for F814w is ~ 7.8, which is larger than the redshift at which 
these peaks are seen. 
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Figure 6. The ^-distribution of 22.5 < bj < 24.0. Panels (a), (b) and (c) are for Scenario A, B and C respectively. Models arc indicated 
by lines. The normalization is the same as Figure 3 to Figure 5. 



3.3.2 z- distribution of 22.5 < bj < 24.0 

Since there are both abundant and deep enough spectro- 
scopic samples in £?-band surveys up to date, the predictions 
of z distributions in this passband is of great significance. 
In Figure 6, we present the predicted z distributions by our 
PLE models within the magnitude range 22.5 < bj < 24.0. 
From panels (a), (b) and (c) we can see that: 1) The high- 
z peaks for Zf = 5.0 models are absent in all of the three 
cosmological models adopted currently. The shapes of these 
curves are similar to those of non-evolutionary ones, except 



that for Scenario A, where there is merely a little rise beyond 
2=3, and 2) the high- z peaks or tails still exist in z/=2.5 
or 3.5 models in any scenario. As mentioned in Section 2, 
the light coming from beyond zl will be greatly reduced by 
Lyman absorption, hence the objects beyond zl are hardly 
detectable, and that zj > 4.0 (for b-band, zl = 4.0) is 
favoured by the z distribution of 22.5 < bj < 24.0 in any 
world models. On the contrary, z/=2.5 or 3.5 are not ap- 
propriate choices, unless the high-z tail or the high-z peak 
can be detected in the future. 
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Figure 7. The differential number counts for E/SO galaxies as a function of apparent magnitude in 7f814 band. The models are 
for Zf = 5.0, which are shown by lines. The letter A, B and C in the figure represent Scenario A, B and C, respectively. Panel (a): 
Comparison with the results of varying the faint-end slope of LF; (b): Comparison with the results of PBZ96. 



4 DISCUSSION 

4.1 Uncertainties in LF 

As mentioned in Introduction, there may exist many uncer- 
tainties in determining the local LF, hence there is no reason 
in believing that the LF in our present work is the unique 
adoption. We will explore the influence of the uncertainties 
in LF upon our conclusions, in particular, upon the num- 
ber counts, simply by varying the LF parameters. From the 
previous investigation, we have already seen that the nor- 
malization, which involves both <j)* and M* (cf. Ellis 1997), 
is appropriate. Thus we only concentrate on the faint-end 
slope a. The following concerns only the z/=5.0 models. 

For Scenario A, we choose a higher value for a, say, - 
0.85, which is slightly lower than -0.9, the value of Marzke 
ct al. (1994a, 1994b). For Scenario B and C, we adopt - 
0.48, the value of PBZ96. From Figure 7-(a) we can see 
that such variation will not affect the counts at bright mag- 
nitudes, and will not bring about too much change at the 
faint-end, either. For Scenario A, such variation leads to 
only 0.12dex higher than that of a = —0.70 at 7fsi4 = 25.0. 
But the faint-end slope of the number-count curve 7 is bet- 
ter than its counterpart of a = —0.70. As for Scenario B 
and C, such variation does not make too much difference. 
Hence, through such prescriptions, we have verified that the 
uncertainties in LF (will) not influence too much our con- 
clusions. 

4.2 Comparison with Other Works 

We notice that the PLE model by PBZ96, which was con- 
structed under the BC93 spectral evolution models in a 



world model with Q ~ 0, can account for most of the ob- 
served photometric and spectroscopic properties of galaxies, 
including the number counts in the U,bj,rf,I and K bands, 
as well as the color and redshift distributions derived from 
most of the existing samples. We have interest in examin- 
ing whether the parameters of the PBZ96 model adopted 
for E/S0 are suitable for the modelling of early- type galaxy 
number counts in the 7fsi4 passband. The e + k- correction 
for E/S0 galaxies in Jfsi4 band is computed by employing 
BC93 models, exactly in terms of the parameters adopted 
by PBZ96. From Figure 7-(b), we can see that, compared 
with our models, neither the t\ nor the T2 model of PBZ96 
(Following PBZ96, n and T2 refer to the time scale of SFR 
T e being lGyr and 2Gyr, respectively, with IMF adopted as 
Scalo one for ellipticals.) can reproduce the number counts 
of E/S0 galaxies properly. Furthermore, the adoption of a 
~ cosmological model is obviously unphysical. 

As can be seen from Section 3.2, our PLE models can 
reproduce the number counts well in all three cosmological 
models under consideration, especially for Scenario C (the 
world model dominated by a A), as seems to be conflict- 
ing with the result of Driver et al. (1996), who concluded 
that flat models dominated by a cosmological constant are 
ruled out from comparison of their E/S0 counts (Driver et al. 
1995b) with their model predictions. The difference between 
our models and Driver et al.'s (1996) lies in the normaliza- 
tion of E/S0 galaxy counts, i.e., a much higher normalization 
adopted by Driver et al. Besides, the E/S0 counts for MDS 
of Driver et al. are apparently lower than the HDF counts 
of Abraham et al. (1996). 
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Figure 8. The differential number counts for E/SO galaxies as a function of apparent magnitude in /f814 band. The models are for 
Zf = 5.0, which are shown by lines. In the two panels, the symbols 'Sc' and 'Sp' denote Scalo and Salpeter IMF, respectively. 'No-ex' 
means no dust extinction is involved. Panel (a): Comparison with the results of no dust extinction; (b): Comparison with the results by 
assuming Salpeter IMF. 



4.3 Dust extinction and IMF 

We have taken into account to some extent the influence 
of dust extinction on the predictions of the models, as de- 
scribed in Section 2.3. From Figure 8-(a), we can see that 
the effect of dust extinction is obvious, especially for Sce- 
nario B and C, such that models without extinction will 
predict more ellipticals than observed from Isi4 ~ 19 to 
~ 21. Moreover, the high-z tails or peaks are even higher 
than the results if dust is not considered in the models (see 
Figure 9). 

The Scalo IMF is more favoured by the observations 
than the Salpeter IMF. It can be seen from Figure 8-(b), 
even with dust extinction involved, the models can not re- 
produce satisfactorily the number counts, and greatly over- 
estimate the counts between I S i4 ~ 19 — 21. A less steep 
IMF at high-mass end such as the Salpeter one will lead 
to more massive stars existing at early times, rendering UV 
fluxes are so strong that more galaxies can be detected at 
high-z (see Figure 9). 

5 SUMMARY AND CONCLUSIONS 

We believe that in any realistic evolutionary models, the 
luminosity evolution must be considered, since there exists 
a well-known fact that galaxies are composed of stars, and 
stars continuously come into birth and evolve into the post- 
main-sequence. The lifetime of stars depends on their mass. 
This leads to the continuous change of the photometric and 
spectroscopic properties of galaxies. Although nE models 



can work sometimes (Driver & Windhorst 1995a), it is ob- 
viously unphysical and can only be treated as baseline for 
comparison with observations and other models. In this 
paper, we have constructed a series of PLE models in the 
three cosmological models to explain the observed number 
counts for elliptical galaxies in MDS and HDF obtained by 
the HST. We summarize our investigations in the following. 

i) Although some of our model colors do not match quite 
well the observed ones, for example, in all cases, the U — B 
color is slightly redder than the observed, the success of the 
BC97 population synthesis model which our present work is 
based upon is still worthy of being affirmed. More recently, 
Steidel et al. (1996) found that high-redshift galaxies are 
also well described by the (BC93) models. It should be 
pointed out that some complexities are neglected in the cur- 
rent prescription of simple PLE models, e.g., we have not 
modelled the evolution of metallicity in our work. There 
exists, however, a so-called age-metallicity degeneracy stat- 
ing that for stellar populations older than 2Gyr, an increase 
in metallicity by a factor of 2 and a decrease in age by a 
factor of 3 results in almost identical optical and near-IR 
colors (Worthey 1994). Therefore, when a formation red- 
shift Zf is assumed so that the age of galaxies is determined 
in a specific cosmological model, the effect of evolution of 
metallicity can be partly compensated by the appropriate 
adoption of the parameters of SFR r e . This will also from 
the other side indicate that the local colors are insensitive 
to the adoption of Zf, as can also be seen from Table 3. 

ii) The Hubble constant Ho plays just a minor role in 
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Figure 9. The z-distribution for ellipticals. Models are indicated by lines. Meanings of the symbols are the same as those in Figure 8. 
Panel (a) is for Scenario A, limited in 22.5 < bj < 24.0; panels (b) and (c) are for Scenario B and C respectively, limited in Ign < 21.0. 



our present investigation. Although Ho is included in the 
normalization cj>* and characteristic luminosity L* of LF, it 
can be canceled out by the quantities relative to cosmol- 
ogy. The only aspect needed to pay attention to is that 
the age of galaxies is measured in Hq 1 . Hence the age of 
galaxies will be smaller if a larger Ho is adopted, and vice 
versa, as will more or less influence the modelled spectra of 
galaxies. The adoption of Ho=50km s _1 Mpc -1 seems quite 
reasonable for a fiat world model with $7 = 1. We can see 
from Table 3 that the integrated colors computed by the lat- 
est BC97 spectral evolutionary models are appropriate com- 



pared with the observed locally (Table 2) . Hence regardless 
of the r7o/globular cluster age problem, the fiat cosmologi- 
cal model will not suffer from the conflict of color/age in the 
context of BC97 models. 

iii) In spite of the fact that the models of Zf = 5.0 failed 
to predict the faint-end flattening of the number counts in 
Scenario B or C, our PLE models can still be regarded to 
have reproduced well the data from the MDS and the HDF 
up to Ifsu ~ 25 within the range of uncertainties in all the 
three scenarios. The corresponding predicted z-distribution 
of 22.5 < bj < 24.0, though we have not observed data to 
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compare with, does not present a high-z tail or peak. The 
shapes of evolutionary predictions are similar to those of no- 
cvolution for the z distribution, with the mean redshift z m 
slightly higher than the no-evolution's, as can be attributed 
to the effect of luminosity evolution. 

iv) The above conclusions are not affected by the variation 
of the LF parameters (0*, a and L*). Hence, even though 
there may exist some uncertainties and biases in the deter- 
mination of local LF, the above conclusions are still valid 
and reliable. 

v) The law of dust extinction for ellipticals is largely un- 
certain. We simply extrapolate the Wang's approach to the 
current case, and the results seem to be well in agreement 
with the observed extinction in B band for present-day el- 
lipticals. By such an ad hoc assumption, together with 
the adoption of Scalo IMF for ellipticals, UV fluxes can 
be greatly reduced at high redshifts, leading to the high- 
z peaks substantially depressed. Moreover, the predictions 
of number counts are also improved. 

vi) High Zf, say 5.0, and nearly-passive evolution shortly 
after a (nearly) single burst to form most of the stars (in 
our work, T e =0.2, and 0.05Gyr for Scenario A, B and C, 
respectively) are critical for the modelling, as is just the 
case of the conventional scenario for formation and evolu- 
tion of ellipticals (cf. Eggen, Lynden-Bell & Sandage 1962; 
Partridge & Peebles 1967). On the contrary, it will not be 
favoured by, especially, spectroscopic samples that galaxies 
formed at as low as Zf = 2.5 (or 3.5), or the evolution is 
as large as that of PBZ96, unless more exotic scenarios are 
introduced into the modelling. 

vii) We do not touch upon the possibilities of number 
evolutions in the present investigation, though it has been 
widely considered to account for the optical/infrared and 
the photometric/spectroscopic paradox (see Introduction). 
Even though there may be the possibility that a new popu- 
lation of dwarf galaxies existed at high-z and then faded or 
disappeared lately, it is unlikely that they are of early-type 
galaxies, or at least most of them can not be early-types. If 
so, the faint-end slope of number counts of early-type galax- 
ies will be steeper than the currently observed data from 
HST. The assumption of mergers is also not necessary here, 
though the possibility of mergers is not ruled out. If the sce- 
nario of merging is taken into account, a series of problems 
will be involved (cf. Toth Ostriker 1992; Babul & Ferguson 
1996, and references therein), and the physical mechanism 
for mergers is much more complicated than that of PLE. 

viii) We have chosen three cosmological models for our 
work. The predictions of number counts by the flat Einstein- 
de Sitter model agree well with the observed data, so do the 
open Friedmann-Robertson- Walker model (f2o=0.1, Ao=0, 
_ffo=50 km s _1 Mpc -1 ) and the Friedmann-Lemaitre model 
(n =0.2, Ao=0.8, #0=60 km s" 1 Mpc -1 ). Therefore, we 
can not discriminate between these cosmological models by 
the galaxy number counts. Our conclusion here contradicts 
with that of Driver et al. (1996), who announced that the 
world models dominated by a A are ruled out by their mod- 
elling of E/SO number counts. Like many other authors 
(Tinsley 1972; Ellis 1997), we believe that the formation 
and evolution of galaxies is a problem far from being well 
understood, and any definite judgment about world models, 
whether the universe is open, flat, closed or A-dominated, 
can not be made as yet by the number counts. 



ix) As mentioned previously, our PLE models do not in- 
clude the evolution of metallicity with respect to z, and those 
details such as the recycling of the residual gas ejected by 
the dying stars as well as some selection effects are also ig- 
nored in our models. These can be considered as further 
improvement over our present models. 

All in all, the models presented in this work can ex- 
plain well the number counts of elliptical and SO galaxies 
derived from HST under the assumption of PLE. It is un- 
likely that the problem of FBGs is caused by early-type 
galaxies, and we hold identical views as other researchers 
(Driver & Windhorst 1995; Driver et al. 1995; Glazebrook 
et al. 1995; Abraham et al. 1996) on this. Furthermore, 
our work shows that the cosmological parameters can not 
be determined by up-to-date observations. In a forthcom- 
ing paper, we will devote ourselves to the galaxy number 
counts of the other morphological types segregated by HST, 
namely, the early-type spirals (Sabc) and late-type spirals 
and irregulars (Sdm/Irr), as well as the overall populations, 
with the present results and conclusions involved, which will 
be examined further. 
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